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ABSTRACT 

Infrared-dark clouds (IRDCs) are believed to be the birthplaces of rich clusters and thus contain the earliest 
phases of high-mass star formation. We use the Green Bank Telescope (GBT) and Very Large Array (VLA) 
maps of ammonia (NH3) in six IRDCs to measure their column density and temperature structure (Paper 1), and 
here, we investigate the kinematic structure and energy content. We find that IRDCs overall display organized 
velocity fields, with only localized disruptions due to embedded star formation. The local effects seen in NH3 
emission are not high velocity outflows but rather moderate (few km s" 1 ) increases in the line width that exhibit 
maxima near or coincident with the mid-infrared emission tracing protostars. These line width enhancements 
could be the result of infall or (hidden in NH3 emission) outflow. Not only is the kinetic energy content 
insufficient to support the IRDCs against collapse, but also the spatial energy distribution is inconsistent with a 
scenario of turbulent cloud support. We conclude that the velocity signatures of the IRDCs in our sample are 
due to active collapse and fragmentation, in some cases augmented by local feedback from stars. 
Subject headings: techniques :interferometric, spectroscopic, stars: formation, ISM: clouds, ISM: kinematics 
and dynamics, radio lines: ISM, Galaxy: structure 



1. INTRODUCTION 

Star formation has been the focus of observational and theo- 
retical studies for decades, but still the conditions under which 
this process commences are quite uncertain. The identifica- 
tion of objects in different evolutionary stages, such that a 
sequence can be constructed, is the essential observational 
ingredient needed to test theoretical scenarios. In the so- 
lar neighborhood, it is possible to resolve the precursors to 
stars (or multiple systems), known as pre-stellar cores, but the 
counterpart in massive regions has to date been difficult to 
isolate. With the recent surveys by Spitzer in the mid-infrared 
and advancement of millimeter and radio interferometric ar- 
rays, progress in identifying objects in various early phases of 
massive star formation has been rapid. 

Infrared-dark clouds (IRDCs), the densest parts of molec- 
ular cloud complexe s embedded within Galactic spiral arms 
( Jack son et al.l |2008[) . are believed to host these earliest 
stage s of clustered star formation Studies in the infrared 
(e.g. Perault et al. 1996; Esanetal. 1998: iRagan et al.ll2009l: 
iButler & T an 2009; Peretto & Fuller 200 91), millimeter cori - 
tinuum (e.g. Ikathborn e et al.l [20 06: Vasvuni na et al.l 12 009). 
and molecular lines (e.g . lCarev e t al. 1998, 2000; Rasan eUri] 
[20061: iPillai et~aTl [20061: ISakai et al.1 12008: Du & Yangll2008h 
have shown that IRDCs contain from tens to thousands of so- 
lar masses of dense (N(H2) ~ 10 22 " 23 cm -2 ) material and have 
the right physical conditions (T < 15 K, n > 10 5 cm -3 ) to give 
rise to rich star clusters, i.e. clusters which can potentially 
host massive (M > 10M o ) stars. 

Star formation is dynamical by nature (see 
McKee & Ostriker 2007, for a review of the important 
processes), but observational tests of dynamics are compli- 
cated by the projection of this three-dimensional process 
onto the two-dimensional plane of the sky. Molecular line 
emission - from a number of molecules excited in the cold 
environments of molecular clouds - is the key tool to help 



disentangle the problem along the line of sight. Ammonia 
(NH3) has been a particul arly useful probe in molecular 
clouds dHo & Townesl Il983|) . as it not only provides kine- 
matic information but also ser v es as a cloud the rmometer 
(IWalmslev & Ungerechtsl 119831: iMaret et al.1 120091) . Am- 
monia has been used widely to s t udy local clouds (e.g. 
iMvers & Benson! 119831; lLadd et al.l Il9g4j; IWiseman & Hoi 
| 1998 | :lJijina etal.l[T 999: Rosolowskv et al. 2008; Friesen et al.1 
2009) and IRDCs fe.g. lPillai et al.ll2006l: iDevine et al.ll20l"iE 
These studies focus on the lower metastable states, (/, K) = 
(1,1) and (2,2), sensitive to the coldest (<20 K) gas without 
any evidence of depleti on. 

In lRagan et al.l (1201 ll hereafter Paper 1), we detailed Very 
Large Array (VLA) observations mapping six IRDCs in the 
NH 3 (/, K) = (1,1) and (2,2). We used the maps to produce 
column density and gas temperature profiles. With between 
4 and 8" angular resolution, we find that ammonia traces 
the absorbing stru cture seen at 8 and 24 [im with Spitzer 
( Rag an etlril 12009). and there is no evidence of depletion of 
ammonia in IRDCs. We estimated a total ammonia abundance 
of 8.1 x 10" 7 and found that the gas temperature is roughly 
constant, between 8 and 13 K, across the clouds. Here, we 
further our analysis of these ammonia data, focusing on the 
velocity structure of the clouds. The high angular resolution 
allows us to profile the kinematics and examine their dynam- 
ical state and stability. 

2. DATA & METHODS 

We obtained observations of the NH3 (1,1) and (2,2) inver- 
sion transitions with the Green Bank Telescope (GBT) and 
Very Large Array (VLA). The observations are described in 
detail in Paper 1. The single-dish and interferometer data 
were combined in MIRIAD (a full description of the method 
is found in Paper 1). A summary of the target properties sen- 
sitivity and resolution of the combined data set is given in 
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TABLE l 

Target summary 



IRDC 




RA 


DEC 


Distance a 


Vlsr 


rms b 


beam size 


Mi RDC c 


Area 


D d 

£>cr 






(J2000) 


(J2000) 


(kpc) 


(km s- 1 ) 


(mJy) 


(" x ") 


(10 3 M ) 


(pc 2 ) 


(mG) 


G005.85- 


-0.23 


17:59:51.4 


-24:01:10 


3.14 


17.2 


2.8 


7.7 x 6.8 


5.5 


0.55 


1.63 


G009.28- 


-0.15 


18:06:50.8 


-21:00:25 


4.48 


41.4 


4.8 


8.3 x 6.4 


1.8 


1.8 


1.43 


G009.86- 


-0.04 


18:07:35.1 


-20:26:09 


2.36 


18.1 


4.3 


8.1 x 6.3 


2.6 


1.3 


0.87 


G023.37- 


-0.29 


18:34:54.1 


-08:38:21 


4.70 


78.5 


2.5 


5.7x3.7 


10.9 


4.4 


1.05 


G024.05- 


-0.22 


18:35:54.4 


-07:59:51 


4.82 


81.4 


4.3 


8.2 x 7.0 


4.0 


1.8 


0.99 


G034.74- 


-0.12 


18:55:09.5 


+01:33:14 


4.86 


79.1 


6.8 


8.1 x 7.0 


5.5 


1.6 


1.43 



a Assume the kinematic near distances from Ragan et al. ( 2006). 

b of the combined data set. 

c from absorption at 8 //m ( Ragan et al. 2009). 
d Critical field strength, see P31 and eq. (0. 

Table \T\ The combined data set has a velocity resolution of 
0.6 km s" 1 . In Table [T] we also list the estimated mass and 
cloud area based o n 8 um extinction, which was computed 
with Spitzer data in Ra gan et al.l (|2009|) and the critical mag- 
netic field stre ngth required for support, which will be dis- 
cussed in 34.31 

At each position, the ammonia spectra were fit with a cus- 
tom gaussian fitting algorithm utilizing the IDL procedure 
gaussf it. The configuration we used for the VLA backend 
did not fit the entire NH3 (1,1) hyperfine signature (spanning 
-3.6 MHz) in the bandpass (3.125 MHz). Our line-fitting 
routine takes a "first guess" line center velocity of the cen- 
tral line (from Rag an et aT1 l2006. see Table [T} which is offset 
by approximately 7.7 km s" 1 from the neighboring hyperfine 
components to either side. We fit each of the components in- 
dependently. For the NH3(2,2) lines, a single gaussian was fit 
to the line independently of the results of the (1,1) fit. From 
these fits, we extract the peak intensity, line-center velocity, 
and Gaussian width of the lines at each position. 

3. RESULTS 

Figures[T]through[6]show the NH3 (1,1) integrated intensity 
and (2,2) integrated intensity plotted over the 8 and 2 4 yum 
Spitzer images of the regions (from iRagan et aI1l2009[) . re- 
spectively, and maps of the line center velocity (first moment) 
and line width (second moment) for the central component 
of the NH3 (1,1) signature. We also plot the physical scale 
assuming the distances in Table [T] 

3.1. Properties of individual sources 

Paper 1 demonstrates that the gas in these IRDCs exhibit 
uniform temperatures, changing by only a few Kelvin in a 
given object, not significantly more than the error. In contrast, 
their velocity fields - both the line-center velocities and line 
width measurements - show a connection between the 
presence of embedded star formation activity and complex 
kinematic signatures. The upper panels of Figures 1-6 show 
the (zeroth, first, and second, from left to right) moment 
maps derived from the NH3(1,1) observations for each IRDC 
in the sample with symbols indicating the locations of the 
24 fim point sources and other young stars, and the lower 
panels show the same for the NH3(2,2) emission. Table [2] 

1 Since our VLA bandpass included only the central 3.125 MHz of the 
~3.5 MHz hyperfine signature, we observe only the central three of the five 
main components of the NH3 (1,1) hyperfine signature. To compute the inte- 
grated intensity of the line, we assume that the missing outer-most lines are 
0.22 the strength of the main lines and that the linewidths of all components 
are equal. 



summarizes the kinematic properties of the emission peaks. 
Although we list only the velocity properties of the central 
component of the NH3(1,1) line, the velocity structure traced 
by the satellite lines closely follows the trends seen in the 
central line. We also list the main line optical depth of 
the NH3(1,1) transition, r m (l,l), the ratio of thermal to 
non-thermal contributions to the pressure (R p ), which will 
be discussed in Section |H and notes about the velocity trend 
in the cloud or particular characteristics of the integrated 
intensity peak. In this section, we discuss the centroid 
velocity and linewidth trends in each IRDC individually and 
connect the detected 24 jim point sources to the kinematic 
signatures. 

For the sake of our modeling, we categorize each IRDC 
based on its NH3(1,1) emission morphology as either a 
"sphere" for objects with an aspect ratio, r, close to one or 
a "filament" for objects with r much greater than one. The 
axes used to make this distinction are indicated in Figures 1- 
6, and the morphological type is listed in Tabled For spheres, 
the aspect ratio is no greater than 1.1, and the elongated struc- 
tures, or "filaments," range from 1.6 to 2.9 in r. 
G005.85-0.23: This source appears approximately round 
(r ~ 1.1) in the NH3 (1,1) and (2,2) integrated intensity map. 
The peak at a(2000) = ll h 59 m 51A\ £(2000) = -24°01'10" 
corresponds to the position of the peak in 8 jim optical depth. 
There are no 24 /jm sources in the mapped region. 

The smooth gradient in centroid velocity in this IRDC 
permits us to straightforwardly quantify and distinguish 
the large-scale ordered motions and the remaining residual 
motion on small scales. We show in the central panels of 
Figure Q] a clear velocity gradient oriented 30 degrees east 
of north. The total gradient in the NH3 (1,1) emission is 
1.2 km s" 1 over 35 arcseconds, or 0.5 pc, resulting in a 
velocity gradient of 2.4 km s" 1 pc -1 . If this linear gradient 
is subtracted, the residual values do no exceed 0.2 km s" 1 , 
indicating the bulk motion dominates the dynamics of the 
cloud. The overall linewidth measured across the cloud is 
very low, between 1.3 and 1.8 km s" 1 , but it increases sharply 
at the edges (to ~3 km s" 1 ) where the centroid velocity also 
falls off quickly. 

G009.28-0.15: In this "filament" (r -1.6), there are three 
integrated intensity maxima: the central peak (a(2000) = 
18*06 m 49.9 5 , £(2000) = -20°59'57", P2 in Table [3, which 
has a 24 [im source associated with it, P3 to the north (off- 
set 25 "), and the maximum (PI) to the south (offset 30 "). 
P2 is near the linewidth maximum (3.3 km s" 1 ), and is also 
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Fig. 1. — Spectral moments in G005.85-0.23. Top left: Spitzer/IRAC 8//m image with NH3(1,1) integrated intensity contours overlaid. Contours begin at 0.2 
Jy beam -1 km s -1 and increase in 0.1 Jy beam -1 km s -1 steps. Bottom left: Spitzer/MIPS 24/um image with NH3(2,2) integrated intensity contours overlaid. 
Contours begin at 0.02 Jy beam -1 km s -1 and increase in 0.01 Jy beam -1 km s -1 steps. Top center: NH3(1,1) centroid velocity map in km s -1 . Top right: FWHM 
of NH3(1,1) central line in km s . Bottom center: NH3(2,2) centroid velocity map in km s -1 . Bottom right: FWHM of NH3(2,2) line in km s -1 . The red line 
represents the major axis and the green line represent the minor axis. The star symbol represents a point source which only appears at 24 jum. The VLA beam is 
shown at the lower-left corner of the first and second moment panels. 



TABLE 2 

Summary of NH3 (1,1) and (2,2) peak characteristics 



Peak J Tdv Position of NH 3 (1,1) transition 



NH 3 (2,2) line 



IRDC 






a 


6 


Vlsr 


Av 


r m (U) 


Vlsr 


Av 


R P 


Notes 


name 






(J2000) 


(J2000) 


(km s -1 ) 


(km s -1 ) 




(km s -1 ) 


(km s -1 ) 






G005.85- 


-0.23 




17:59:51.4 


-24:01:10 


17.36+0.02 


1.41+0.03 


4.9 


17.49+0.04 


1.37+0.04 


0.1 


smooth v-grad. 


G009.28- 


-0.15 


PI 


18:06:50.8 


-21:00:25 


40.99+0.02 


1.45+0.02 


4.0 


40.99+0.03 


1.53+0.03 


0.1 


main peak 






P2 


18:06:49.9 


-20:59:57 


41.70+0.07 


2.37+0.07 


4.0 


41.75+0.08 


2.14+0.08 


0.04 


24 fim source 






P3 


18:06:49.8 


-20:59:34 


41.47+0.03 


1.85+0.03 


3.6 


41.38+0.04 


1.67+0.05 


0.06 


24 yi/m source 


G009.86- 


-0.04 




18:07:35.1 


-20:26:09 


17.75+0.04 


1.35+0.04 


3.3 


17.50+0.06 


1.64+0.07 


0.1 


"quiescent" peak, 2 v-grad. 


G023.37- 


-0.29 




18:34:54.1 


-08:38:21 


78.81+0.20 


3.89+0.22 


a 


78.16+0.15 


3.66+0.19 




24 11m source 


G024.05- 


-0.22 




18:35:54.4 


-07:59:51 


81.65+0.03 


1.96+0.03 


2.6 


81.62+0.05 


2.14+0.07 


0.05 


N-S v-grad. 


G034.74- 


-0.12 


PI 


18:55:09.5 


+01:33:14 


77.95+0.03 


2.46+0.03 


6.1 


77.72+0.05 


2.08+0.05 


0.04 


24 11m source 






P2 


18:55:11.0 


+01:33:02 


78.69+0.03 


2.09+0.03 


2.9 


78.58+0.07 


2.04+0.07 


0.05 


24 jim source 


a line saturated. 



red-shifted in centroid velocity. PI, while the strongest in in- 
tegrated intensity has the lowest linewidths detected in this 
object (1.4 km s" 1 ), and there is no associated 24 /mi source. 
The northern integrated intensity peak is 10 " away from a 
24 yum point source, but the kinematic structure is not altered 
by its presence. 

Apart from P2 and P3, the bulk of the cloud resides at a 
narrow range of line center velocities, between 41 and 41.5 
km s" 1 . The sharpest changes in centroid velocity are located 



at the eastern edge of the cloud, where the line is blue- shifted 
by 1-1.5 km s" 1 with respect to the bulk of the cloud at the 
southeast edge in both NH 3 (1,1) and (2,2) emission. The 
linewidths are also enhanced at this edge, though no YSOs 
are detected in this region. 

G009. 86-0.04: We approximate this source as a filament, the 
most elongated structure (r ~2.9) in our sample. The inte- 
grated intensity peak (a(2000) = 18 /z 07 m 35.0', £(2000) = 
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Fig. 2.— Same as Figure[]]but for G009.28-0.15. 

-20°26'09") is dark at both 8 and 24 fim and corresponds 
to where the centroid velocity and linewidth is the lowest, all 
evidence for a quiescent region. The NH3 (1,1) centroid ve- 
locity field in this object is organized into two gradients in 
either direction from the central peak in integrated intensity. 
The eastern (left-hand) gradient is of magnitude -1 km s" 1 
oriented 80° east of north, and the western (right-hand) gra- 
dient is of magnitude -1.5 km s" 1 oriented 65° west of north. 
The central "hinge" position is indistinct in the linewidth mea- 
surement. 

Overall the velocity field in the filament is smooth, with 
several YSOs coincident with the ammonia emission: five 
east of the intensity peak, and one to the southwest. The 
NH3 (1,1) linewidths are enhanced (> 2 km s" 1 ) in the east. 
Curiously, the NH3 (2,2) emission, which appears to follow 
the locations of the YSOs in the eastern region, exhibits an 
overall shift to lower line-center velocities (by 0.6 km s" 1 ) 
and higher linewidths (by 0.5 km s" 1 ). The optical depth of 
the NH3(1,1) main line is below -3 in the eastern region, so 
it is unlikely that optical depth effects are the cause of the 
increased linewidth. This object appears to be undergoing 
cluster formation, though the part of the cloud associated 
with NH3 (1,1) emission peak remains quiescent. 

G023.37-0.29: The integrated intensity map of this 
round IRDC (r -1.1) peaks at a(2000) = 18*34 m 54.1 5 , 
£(2000) = -8°38'21", although throughout this cloud, 
the lines are saturated and/or optically thick, making it 
impossible to derive reliable optical depths or very accurate 
line properties. There is a 24 jim point source (not present at 
8 jjm) near the center of the region, slightly offset from the 
intensity peak, which is likely a deeply embedded protostar. 



The central region has very high linewidths (highest in the 
sample, 4 km s" 1 , see Figure |4]), although because of the 
high optical depth of the NH 3 lines, these should be taken 
cautiously. 

G024.05-0.22: This approximately round source (r -1.1) 
appears centrally-peaked in line intensity (left panels of 
Figure [3 at a(2000) = 18 /z 35 m 54.P, £(2000) = -7°59'51", 
which corresponds also to the peak in 8 fim optical depth. 
The NH3 (1,1) maps shows a velocity gradient starting at an 
east- west aligned "ridge" slightly offset to the north from the 
peak of integrated intensity. This "ridge" in centroid velocity 
also corresponds with enhanced linewidths (-3.6 km s" 1 ) in 
both the (1,1) and (2,2) lines, though with no distinction in 
integrated intensity similar to what we see in G009. 86-0.04. 
From the center of the cloud across this ridge, the velocity 
changes by -1.1 km s" 1 , and corresponds to a gradient of 
2.1 km s" 1 pc" 1 . At the southern tip of (1,1) emission, there 
appears to be a clump with distinct red- shifted velocity but 
this is not detected in (2,2) emission. While there is one 
Class II source coincident with northern ridge, this IRDC 
lacks 24 jjm detections (an indicator of an embedded source) 
anywhere in the cloud. 

G034. 74-0.12: The overall velocity field of this filament 
(r -2.0) appears quite disorganized, particularly in locations 
where 24 jum point sources are detected. There are two in- 
tegrated intensity peaks in this IRDC, both in the vicinity 
of 24 yum sources (star symbols in Figure 0), and both po- 
sitions show high linewidth. The strong NH3 (1,1) peak in 
the northwest portion of the cloud (a(2000) = 18*55 m 09.5 5 , 
£(2000) = + 1°33'14", PI in Table [2]) is directly coinci- 




Fig. 3. — Same as Figure[T]but for G009. 86-0.04. NH3(2,2) integrated intensity contours (overplotted in bottom-left panel) begin at 0.03 Jy beam 1 km s 1 and 
increase in 0.01 Jy beam -1 km s -1 steps. 



dent with a 24 //m point source, enhanced linewidths and a 
slightly blue- shifted centroid velocity. The optical depth of 
the NH3(1,1) main line is 6.1. The centrally located peak 
(a(2000) = 18^55^11.0^, £(2000) = +1°33'02", P2 in Ta- 
ble [2]) is offset 10 " from the position of the 24 jum source 
and offset 15 " from the nearby peak in linewidth. The optical 
depth of the NH3(1,1) line here is 2.9. These two locations 
are also where most of the appreciable NH3 (2,2) emission is 
found. 

3.2. Comparison between (1,1) and (2,2) kinematics 

As is shown in the left panels of Figures [T] through [6] the 
NH3 (1,1) emission tends to be more widespread than the (2,2) 
emission. In this section, we compare the velocity fields of the 
two states. The central panels show the range in line center ve- 
locities, which generally encompass the same range, and the 
right panels show similar trends in linewidth. Table [2] shows 
the velocity properties of the (1,1) and (2,2) at the locations 
of the NH3(1,1) intensity peaks. 

In G009.86-0.04 and, to a lesser extent in G023.37-0.29, 
we find that the median NH3(2,2) line center velocity is offset 
to lower velocities by -0.6 km s" 1 compared to the NH3(1,1), 
and the median FWHM of the (2,2) line is higher than that of 
the (1,1) line by ~0.8 km s" 1 . Figure [7] shows the distribu- 
tions of line center velocity and linewidth from both the (1,1) 
and (2,2) maps. G009. 86-0.04 has several young stars and 
24 yum point sources in the eastern part of the cloud, which 
corresponds to the locations of the high linewidths and blue- 
shifted line center velocities in the NH3(2,2) map. The optical 
depth of the NH3(1,1) main line is less than 3 in this region, 
lower than typical values of 4 or 5 throughout the sample, so 



it is unlikely that optical depth effects are the cause of the en- 
hanced linewidth. It may be the case that the young stars are 
having a dynamical effect on the slightly warmer gas probed 
by the NH3(2,2) emission. This does not appear to be the case 
in G023.37-0.29, where there is only a singular 24 jim point 
source and the lines are optically thick, which would likely 
limit our ability to probe near any embedded source(s). 

IRDC G034.74-0.12 also exhibits evidence that a cluster 
is actively forming with the two 24 jim point sources cor- 
responding to very strong NH3 emission and several other 
young stars in the vicinity. However, in these two locations 
of 24 yum point sources, which is also where most of the 
NH3(2,2) emission is detected, the linewidth appears more 
enhanced in the NH3(1,1) rather than (2,2). We note that 
the high optical depth of NH3(1,1) at PI (6.1) may contribute 
to the broadened line here. The observations of this object 
were the noisiest of the sample, and it is also the most dis- 
tant IRDC, so we may not be sensitive to the effect we see in 
G009.86-0.04 (the nearest IRDC). 

3.3. Linewidth 

The NH3 (1,1) linewidths in our sample of IRDCs are 
between 1.1 and 4 km s" 1 , occupyi ng the high tail o f the 
linewidth distribution presented in the Jiji na et aO |l999) sur- 
vey of 264 de nse cores, but on p ar with other ammonia stud- 
ies of IRDCs: Pill ai et all (|2006), wh o found a s l ightly lower 
range in their single-dish study, and lWang et al.l (l2008|) . who 
found that different cores within an IRDC exhibited differ- 
ent linewidths in the high-resolution study: higher linewidths 
near locations of embedded star formation activity and lower 
linewidths in quiescent regions of the cloud. For our sources, 
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Fig. 4.— Same as Figure[]]but for G023.37-0.29. 

the enhanced linewidth appears to correspond to the locations 
of 24 jim point sources. 

We find that the linewidth increa ses with increasin g dis- 
tance to IRDCs, as was noted by iPillai et al.l (|2006). In 
Figure [8] we show that the maximum linewidth detected in 
IRDCs varies directly with distance, which may be a result 
of clumping within the beam increasing with distance. Cer- 
tainly, with Spitzer we do see objects on the 2 - 3 " scale which 
would not be resolved with the beam (sometimes 6 - 8"in low- 
elevation sources). At the same time, the minimum linewidth 
does not show any trend with distance. As noted above, the 
integrated intensity peaks with 24 jim point sources have a 
higher linewidth than those peaks without, but the apparently 
starless peaks typically have very low linewidths, independent 
of distance (Av -1.4 km s" 1 ). Therefore, in the following 
discussion, we proceed with our analysis assuming that the 
trends are intrinsic features not determined by the distance. 
As shown in the lower panel of Figure [SJ the cloud masses in- 
crease slightly with distance, which could simply be a selec- 
tion effect. Higher cloud masses suggest deeper gravitational 
potentials, which in turn would give rise to larger linewidths. 
Yet the total kinetic energy calculated from the averaged line 
widths and the centroid velocities does not change perceptibly 
with distance. 

The dynamical properties of all IRDCs are summarized in 
Figure [9] The centroid velocity distributions (left column) are 
asymmetric, with tails to lower or higher velocities than the 
systemic velocity. This could indicate a substantial elongation 
along the line-of- sight, or highly asymmetric infall. The rel- 
ative motions to the systemic velocity are mostly supersonic. 



The non-thermal velocity dispersion 

[a^ ^7 

°~ NT = a o~n^> 0) 

V 8 In 2 \mih 

is larger than the sound speed (see Paper 1 for temperatures, 
ranging from 8 to 13 K) for all IRDCs - in fact, none of 
the velocity dispersions are smaller than Mach 2. The dis- 
tributions are peaked, with tails to high Mach numbers. For 
IRDCs containing sources (see discussion above), these tails, 
or "excesses" are correlated with locations in the vicinity of 
the sources. 

4. DISCUSSION 

Judging from Figures [T] through [6l it is unlikely that one 
model of IRDC kinematics can fully account for the range 
of behaviors observed. The broad characteristics - the over- 
all linewidths and ranges of centroid velocities, see Table [2] 
- are roughly consistent among the clouds in the sample and 
with previous mo lecular line studies (e.g. IPillai et al J 12006: 
Ragan et al. 2006). However, our sample exhibits both global 
trends (e.g. smooth velocity gradients) and localized effects 
(e.g. signatures of feedback from young embedded proto- 
stars) that can now be investigated with high angular reso- 
lution observations. 

As a first approximation of stability, we compute the virial 
mass of the clouds (My iria i = 5RV^ ms /(3G)), where R is the 
radius of the cloud, G is the gravitational constant, and V rms 
= 3 2 AV/2.35 where AV is the average linewidth of the cloud. 
Virial masses for the whole clouds are typically 10 2_3 M o , and 
the virial parameters, a = Mviriai/M, from 0.1 to 0.7, suggest- 
ing that IRDCs are bound structures prone to collapse. The 
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Fig. 5. — Same as Figure[T]but for G024.05-0.22. NH3(2,2) integrated intensity contours (overplotted in bottom-left panel) begin at 0.04 Jy beam 1 km s 1 and 
increase in 0.01 Jy beam -1 km s -1 steps. 



cloud masses M ar e taken from the dust extinction maps by 
lRaganetal.ld2QQ9h . 

The linewidths observed in our sample (1.1- 4.0 km s" 1 ) are 
in excess of the thermal linewidth, -0.18 km s" 1 . The ratio of 
thermal to non-th ermal pressure in the cloud, as expressed by 
lLada et ail (|2003|) . is R p = c 2 s lcr 2 NT , where c s is the isothermal 
sound speed, and cr NT is the three-dimensional non-thermal 
velocity dispersion. We calculate values for R p at each of the 
intensity peaks in the IRDCs (see Table [2]) and find an average 
value of 0.07, indicating that non-thermal pressure is dom- 
inant. "Non-thermal effects" encompass many things, such 
as infall, outflows, or systematic cloud motion s (i.e. rotation) 
and p ossible "support" fro m turbulent motions ( Arons & Maxl 
119751) or magnetic fields (iMouschovias & Spi tzer 1976]). In 
the following sections, we explore these effects individually, 
so to determine the dominant processes. 

4. 1 . Connecting star formation and IRDC kinematics 

The IRDCs in our sample span a range of stages of star for- 
mation, some devoid of embedded sources and some hosting 
several. In Paper 1, we showed that the presence of young 
stars does not significantly affect the kinetic temperature of 
the gas traced by ammonia, at least beyond the ~ 1 K errors 
we estimate^. However, as we showed in Section 3. 1 , the pres- 
ence of young stars, particularly 24 /im point sources, has a 
localized effect on the IRDC dynamics, namely an increase in 

2 We noted in Paper 1 that these observations of the (1,1) and (2,2) tran- 
sitions of ammonia are not necessarily sensitive to hotter gas on the small- 
est scales, which (for example) may arise from a young embedded protostar 
heating a compact core. Gas in warm, compact regions are better probed with 
higher- J transitions of ammonia or other molecules. 



linewidth at the positions of young stars and (sometimes) dis- 
tinct centroid velocity components. We see no strong trend for 
the positions of the 24 [im point sources to have exceptionally 
high optical depths in NH3(1,1), thus we do not expect opti- 
cal depth effects to contribute strongly to linewidth enhance- 
ments. Furthermore, starless peaks in NH3 intensity (existing 
in all clouds except G023.37-0.29 an G034.74-0.12) tend to 
have the lowest linewidths. 

Broadened linewidths are an indicator of increased inter- 
nal motions which accompany the onset of star formation 
(iBeuther et al.ll2005l) . In each IRDC (except G005.85-0.23 
and G024.05-0.22) of our sample, we see that the sites of 
24 /mi point sources are accompanied by an enhancement 
in NH 3 (1,1) linewidth (in G009. 86-0.0 4, the enhancement is 
rather seen in the NH 3 (2,2) linewidth). Beuth er et al.l £2005) 
found that in a molecular core with an infrared counterpart, 
the N2H + (l -0) emission (known t o trace dense gas similarly 
to NH3, e.g. Johns tone et al.ll20Tol) exhibited broad line emis- 
sion, whereas linewidths in infrared-dark cores - presumably 
at an earlier evolutionary stage - were significantly narrower 
(Av ~ 1 km s" 1 ). The increased internal motions giving rise 
to broadened lines can be attributed to ordered motion, such 
as infall, outflow or rotation. 

4.2. Kinematics of "starless" IRDCs 

Because the feedback from embedded protostars can con- 
fuse the kinematic signatures in the IRDCs, here we take 
a closer look at the IRDCs which show little or no evi- 
dence of actively forming stars. Our sample includes one 
IRDC, G005. 85-0.23, which lacks any coincident Spitzer 
point sources, thus does not appear to host embedded star 
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Fig. 6. — Same as Figure[T]but for G034.74-0.12. NH3(2,2) integrated intensity contours (overplotted in bottom-left panel) begin at 0.04 Jy beam 1 km s 1 and 
increase in 0.01 Jy beam -1 km s -1 steps. 

TABLE 3 

Summary of fit results 



IRDC 




geometry a 


n b 


To c 


en 

bn 


Q~res 


A f 


name 






(cm -3 ) 


(K) 




(km s -1 ) 




G005.85- 


-0.23 


sphere 


1.1 x 10 6 


5.0 x 10 3 


1.1 


2.4 


0.18 


G009.28- 


-0.15a 


filament 


4.0 x 10 5 


5.4 x 10 3 


2.3 


2.5 


0.23 


G009.28- 


-0.15b 


filament 


8.9 x 10 5 


1.2 x 10 4 


5.5 


3.8 


0.27 


G009.86- 


-0.04 


filament 


6.8 x 10 5 


2.0 x 10 3 


6.2 


1.6 


0.28 


G023.37- 


-0.29 


sphere 


6.0 x 10 5 


6.0 x 10 3 


1.7 


2.6 


0.31 


G024.05- 


-0.22 


sphere 


5.5 x 10 5 


6.0 x 10 3 


1.7 


2.7 


0.4 


G034.74- 


-0.12 


filament 


3.4 x 10 5 


2.3 x 10 4 


2.0 


5.3 


0.18 



a Fitting geometry (BE-sphere or isothermal cylinder). 
b Fitted central density. 
c Fitted isothermal temperature. 

d Normalized stability parameter. % n = £ m \ for cylinders, and^ n = ^5^/6.5 for BE-spheres. % n > 1 indicates instability. 
e Residual velocity dispersion needed at minimum to support the cloud energetically. 
f Normalized rms error of fit. 



formation. G024.05-0.22, has one Class II object (with no 
24 /mi counterpart, see Figure near the edge of the NH3 
(1,1) emitting region, thus the bulk of the IRDC appears de- 
void of YSOs. It is possible in both cases that embedded pro- 
tostars in the IRDC are heavily extincted by dust beyond our 
detection limit, but we continue our discussion assuming that 
the kinematics are dominated by the global forces rather than 
protostellar feedback. 

Both of these IRDCs have nearly round projected mor- 
phologies. In G005.85-0.23, the gradient is from southwest 



to northeast (33° east of north) centered on the integrated in- 
tensity peak. In G024.05-0.22 the gradient is from northeast 
to southwest (8° east or north) but is not symmetric about the 
peak in integrated intensity. In this case, the velocity gradient 
is not smooth, but has a sharp ridge structure in the east- west 
direction, which also corresponds with enhanced linewidths 
(~3.6 km s" 1 compared to < 2 km s" 1 throughout the rest of 
the cloud). 

Smooth velocity gradient s are often interpre ted as 

signatures or rotation (e.g. Arquilla & Goldsmith 1986; 
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Fig. 7.— Distribution of velocity fit parameters for IRDC G009.86-0.04 
and G023. 37-0.29 with the NH 3 (1,1) plotted in the black histogram and the 
NH3(2,2) distribution plotted in the red histogram. Both histograms are nor- 
malized to the number of positions with a (1,1) measurement. The left panel 
shows the line center velocities and the bottom panel shows the FWHM, both 
in km s _1 . 



iGoodman et al.lfl993|) . Since the projected geometry of both 
of these IRDCs is roughly circular, we can reasonably approx- 
imate the clouds as spheres. If we assume solid-body rotation, 
the resulting velocity gradients are 2.4 and 2.1 km s" 1 pc" 1 for 
G005. 85-0.23 and G024.05-0.22, respectively. If this orga- 
nized motion is linearly fit and subtracted from the centroid 
velocity field, the residuals are less than 0.2 km s" 1 . 

We can then compare the importance of rotational kinetic 
energy to gravitational energy, parameterized by the ratio, 
J3 (in this case for a uniform density sphere), defined as 
J3 = £l 2 R 3 /(3GM), where D is the angul ar velocity, R is th e 
cloud radius, and M is the cloud mass ( Good man et al.[[T993h . 
For the mass, we take values from Raga n et alJ (120091) using 
8 /mi absorption as a mass-tracer over the region mapped in 
ammonia. For reference, a value of /3 = | is equivalent to 
breakup speed for a spherical cloud, and lower values signify 
a lessening role of rotation in cloud energetics. Under these 
assumptions, we find values of 2 x 10" 4 and 5 x 10" 4 , lower 
than the extremely low end of the j3 range seen in dark cores 
( Goodma n et al.l fl993|) because the masses are much higher 
(500 and 2500 M© respectively). Adopting a more realistic 
centrally peaked density profile, p oc r" 2 for example, reduces 
J3 by a factor of 3. In this simplistic picture, rotation plays 
only a small role in the dynamics of the cloud. 

Such a simple solid-body rotation model for such spheri- 
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Distance (kpc) 

Fig. 8. — Top: Linewidth vs. IRDC distance. The linewidth at the in- 
tegrated intensity peaks (marked with asterisks where there 24 //m point 
sources present and + signs for those without 24 //m point sources) and the 
maximum detected linewidth overall in the cloud (marked with diamonds) 
both incr ease with inc reasing distance. Bottom: Total mass (from 8 jim ab- 
sorption, Ragan et al. 2009) vs. distance. 

cal IRDCs is a temptin g interpretation, but one that should 
be made cautiously. Burkert & Bodenheimer (2000) show 
that large-scale turbulent motions (i.e. turbulence with a 
steep power spectrum) can lead to centroid velocity gradi- 
ents that look like shear or rotation. Indeed, even in these 
IRDCs which appear to be the most quiescent, the linewidths 
far exceed the thermal sound speed in the typical IRDC en- 
vironment (about 0.2 km s" 1 ). Such linewidths could be 
caused by outfl ows from a low-m ass stellar component not de- 
tectable in the Ragan et al. ( 2009) Spitzer observations. Even 
in the absence of stellar feedback, non-thermal linewidths 
woul d be expected as consequence of the cloud formation 
fe.g. lVazquez-Semadeni et al.ll2007l: iHeitsch et al.ll2008bl) . or 
if glo bal gravitational collapse dominates the clo ud evolution 
(e.g. lBurkert & Hartmannl2004l:lField et al.ll2QQ8h . Below, we 
examine these possibilities in greater detail. 

4.3. Dynamical conditions of IRDCs 

To refine our energetics estimates of the previous sec- 
tion, we analyzed the spatial energy distribution within the 
clouds. To date, dynamical studies of IRDCs have been lim- 
ited mainly to single-dish surveys with resolution elements 
of -30", which is insufficient to resolve the relevant (sub- 
parsec) scales. With our VLA dataset, we have mapped the 
velocity field across an IRDCs, so we are now able to quan- 
tify the energy distribution in IRDCs. To determine the degree 
of stability, we will fit (idealized) geometries, guided by our 
classification of the IRDCs in "spheres" and "filaments" (see 
Tab. [3]). The fits result in density profiles, isothermal temper- 
atures, and a criticality parameter (see below). We use the 
fitted temperatures as a measure of the energy required to bal- 
ance gravitational and kinetic energy content, which can be 
compared to the measured temperatures from Paper 1 . The fit 
results are summarized in Table [3] 

4.3.1. "Filaments" 
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Fig. 9. — Summary of IRDC kinematics. Left column: Histogram of the centroid velocity in units of the sound speed. Vertical dashed lines indicate the transition 
between sub- and super-sonic. Right column: Histograms of the non-thermal velocity dispersion (eq.[T}- All values are supersonic by at least a factor of 2. 



For the filament-like clouds, we construct radial mass den- 
sity profiles from the dust column densities assuming an 
isothermal cylinder as the underlying model. The density pro- 
file for an isothermal cylinder is given by 



with the cylinder scale height 



H 2 = 



7rGp ' 



(3) 



p(r) = 



Po 



(1 + (rlHf)- 



2\2' 



(2) 



(e.g. lOstrikerfl 964) with the sound speed c s = y/kBTo/Qimn). 
The fitting is done in two steps. First, we construct a filament 
following th e mass distribution as traced by dust extinction 
with Spitzer ( Raga netaTI 120 09) by calculating the center-of- 
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mass positions along the RA and DEC axes. A linear re- 
gression through the resulting positions results in the filament 
axis, from which the distances of the sample positions are cal- 
culated. This gives us a radial column density profile N t s (r). 
Next, for the fitting, we start with an initial guess of po, To 
and Ro, construct a profile using equation © and project it to 
generate a column density profile N C yi( r )- The rms difference 

Arms = Yj(N obs (ri) - Ncyiin)) 2 (4) 

i 

between the data and this profile is then minimized by a down- 
hill simplex method in the three parameters po, To, and Ro. 
For each cloud, we check a map of A rms to ensure that the fit 
does not converge on a local minimum (it never did). 

To estimate the degree of stability, we compare the fit- 
ted temperature, To, which range between 2.0 xlO 3 and 2.3 
xlO 4 K, to the cloud temperatures presented in Paper 1 - typ- 
ically between 8 and 13 K. We convert the difference in tem- 
peratures to a "residual velocity dispersion" (cr res ) that would 
be needed to support the cloud against collapse, from 1.6 to 
5.3 km s" 1 . These parameters are summarized in Tabled We 
also take the ratio of the mass per unit len gth derived fro m the 
fit over the corresponding critical value ( Ostriker 1964), 

_ m _ nGRlpo 16H 2 
Ul = m c ~ 2c] 16^+^' (5) 

with a mean molecular weight of \i - 2.36. Note that the fitted 
temperature To is not the temperature needed to stabilize the 
filament against collapse, since we fit all three parameters, 
central density, temperature, and radius. Thus, for £ m / > 1, 
the filament will collapse. 

Figure ITOl summarizes the results. The left column shows 
the dust column profiles (symbols) and the best fit (solid line) 
with the corresponding parameters. As is already suggested 
by the maps, there is a substantial scatter in the profiles, re- 
sulting in relative rms errors between 20 and 30%. Note that 
we split up G009.28-0.15 into two components divided at 
6 = -21° 00' 20", just as was done in Paper 1. The non- 
thermal velocity dispersion profiles (Equation [U are given in 
the right column, with the data shown in symbols, a linear re- 
gression in a dashed line, and a binned version in the solid line 
(including error bars). The dot-dashed line indicates the resid- 
ual velocity dispersion that would be (at minimum) required 
for cloud turbulent support given the nominal fit temperature 
and the (observed) actual temperature. 

All "filaments" are gravitationally unstable by at least a fac- 
tor of 2 in terms of masses, consistent with the virial estimates 
made above. The non-thermal velocity dispersions cr NT , typ- 
ically between 0.5 and 2 km s" 1 (corresponding to hundreds 
of Kelvin), are at least a factor of 2 below the formally re- 
quired value for support, cr res . If cr NT were to give rise to an 
isotropic pressure as envisag ed in models of turbulent support 
(Mc Kee& r j^nll2002L 120031) . then its gradient is inconsistent 
with such a scenario for at least two IRDCs (G023. 37-0.29 
and G009. 28-0. 15b). The remaining IRDCs show strong 
scatter in the velocity profiles, indicating that there is no sys- 
tematic velocity distribution. 

4.3.2. "Spheres" 

IRDCs with roughly circular shapes (or at least not ob- 
viously filamentary ones) we approximate by a ( project ed) 
Bonnor-Ebert sphere. The same reasoning as in 34.3.11 ap- 



plies, namely that we are interested in the nominal temper- 
ature required to provide pressure support agai nst collapse . 
We solve the modified Lane-Emden equation (|Ebertl [19551: 
lBonnorf[l956h 

% = — ^AnGnoumH, (7) 

c s 

integrating the ordinary differential equations with a 4th or- 
der Runge-Kutta scheme, with initial conditions if/(0) = 
and dif/(0)/dt; = 0. A solution is determined by the central 
core density no, the isothermal temperature To, and the radius 
Rmax- For R > R max , the density is assumed to drop to irrele- 
vant values, with the temperature increasing to provide pres- 
sure equilibrium at R max - For values of no/n(R max ) > 14.3, or 
%be > 6.5, the BE-sphere is gravitationally unstable and will 
collapse. 

In a procedure similar to the filament- fitting, we construct 
radial column density profiles from the dust data, which we 
then fit with a two-dimensional projection of a BE-sphere. We 
chose to fit all three parameters, no, To and R max , instead of 
constraining the fits by the observed cloud radius R b s and the 
central column density, because having the BE-sphere extend 
to R max = Robs results in compact column density profiles, 
i.e. the column density would drop to zero at R b s , inconsis- 
tent with the shape of the observed column density profiles. 
Thus (again), the resulting BE temperatures are not the tem- 
peratures needed to support the sphere, but they are generally 
smaller. 

Figure [TT1 summarizes the results analogously to Figure [T0l 
We show the radial column density profiles, the best fit and its 
parameters, and the velocity profiles with a linear regression, 
a binned profile, and the turbulent residual velocity needed to 
support the cloud. Comparing this to the actual dispersions 
gives us a measure of cloud stability. Again the fitted temper- 
atures are on the order of 5-6 xlO 3 K. 

All fits result in unstable BE-spheres. The ^^-values are 
super-critical even at the nominal (fitted) temperatures which 
are more than two orders of magnitude larger than the ob- 
served temperatures. The cr NT -profiles are more than a factor 
of 2 below the dispersion required for support, and the ob- 
served values show again a substantial scatter, indicating a 
strong degree of anisotropy in the spatial kinetic energy dis- 
tribution. 

4.3.3. Implications of fitting results 

The above discussion leads us to conclude the following. 
All IRDCs in our sample are gravitationally unstable, with 
gmi > 1, %be > 6.5 (see Tab. [3]). The fitted temperatures are at 
least two orders of magnitude higher than the observed values 
for all IRDCs, yet the fits still result in unstable structures. 
Interpreting these temperatures as "turbulent temperatures", 
the kinetic energy (including systematic motions as indicated 
by the centroid velocity) in the clouds is insufficient to pro- 
vide support against collapse, consistent with our earlier virial 
estimates. The scatter of the velocity profiles is substantial 
(in cases more than 100%), suggesting that local gravitational 
motions dominate over "micro-turbulent" motions. The latter 
have been repeatedly demonstrated to b e inconsistent with ob- 
servational and numerical evid ence (e .g. lBrunt & Heyeri2002t 
iPadoan et al.ll2003l: iBrunt et alJl2009h . 
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Fig. 10. — Left: Observed column density profiles (symbols) and best fit (line) for the three "filamentary" IRDCs, with G009.28-0.15 split up in a North (a) and 
South (b) component. Central density no, isothermal temperature To and the ratio of mass per unit length over critical mass £ are indicated for each fit. Right: 
Corresponding non-thermal velocity dispersion profiles for each IRDC (symbols). A linear regression is indicated by the dashed line, and a binned profile by the 
solid line including error bars. The dot-dashed constant line shows the "turbulent residual velocity" necessary to stabilize the cylinder. From the values and the 
difference between actual and residual dispersion we see immediately that all IRDCs are unstable. 
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Fig. 11. — Left: Observed column density profiles (symbols) and best fit (line) for the three "spherical" IRDCs modeled as Bonner-Ebert spheres. Central density 
no, isothermal temperature To, and %be are indicated for each fit. Right: Corresponding non-thermal velocity dispersion profiles for each IRDC (symbols). A 
linear regression is indicated by the dashed line, and a binned profile by the solid line including error bars. The dot-dashed constant line shows the "turbulent 
residual velocity" necessary to stabilize the BE-sphere. From the values and the difference between actual and residual dispersion we see that all IRDCs are 
unstable. 




Fig. 12. — Normalized centroid velocity histograms of collapsing and fragmenting molecular clouds (see text). The measurements are centered on the three 
most massive gravitationally bound cores (numbers 2, 5 and 8 from model Gf2 of Heitsch et al. 2008a). Asymmetries and high-velocity tails are visible similar 
to the observational histograms. While not proof, the consistency of distributions is suggestive. 
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Another way to view our results is to consider the "turbu- 
lent residual velocity dispersion" (dot-dashed line in Figs.fTOl 
and [TT]) as a measure for the gravitational energy within 
the IRDC. Then, its ratio to the observed dispersion values 
of a few is consistent with the ratios between gravitational 
and kin etic energies found in models of collapsing molecula r 
clouds ( Vazquez-Semadeniet al. 20071: iHeitsch et a l. 2008a), 
in which the non-thermal ("turbulent") line widths in molecu- 
lar clouds are driven by global gravitational collapse, with the 
kinetic energy trail ing the gravitational energy by a factor of a 
fe w (see Fig. 8 o f IVazquez- Semadeni et al. 2007 and Fig. 10 
of IHeitsch et"aT1l2008al) . 

Expanding on this thought, Figure[T2]shows centroid veloc- 
ity histograms taken along three lines-of- sight from a model 
of flow-driven cloud formation (model Gf2 of Heitsch et al. 
2008a). The centroid velocities were measured at a point 
when the cloud is gravitationally collapsing, and forming lo- 
cal gravitationally bound cores. The selected lines-of- sight 
are centered on three of the most massive cores in the Heitsch 
et al. study (numbers 2, 5, and 8). The clouds form due to the 
collision of two warm, diffuse gas flows. Strong hydrodynam- 
ical and thermal instabilities lead to immediate fragmentation 
and, once a sufficiently high column density has been assem- 
bled, to local and global gravitational collapse. The model 
centroid profiles show similar asymmetries and tails as we 
see in Figure [9] These asymmetries arise from infall in a 
non-uniform medium, i.e. clumps of gas are falling into the 
gravitational potential well. Such events may cause the asym- 
metries in the observed profiles. In other words, gravitational 
collapse of a molecular cloud would not necessarily result 
in symmetric centroid velocity distributions. Further high- 
resolution studies, using appropriate high-density gas tracers, 
must be conducted in order to test our inferences about gravi- 
tational collapse in more detail. 

If turbulence cannot support the IRDCs, could magnetic 
fields? In the absence of observational data for our IRDCs, 
we can estimate the critical field strengths required for cloud 
support. Using the expressi on for the critical mass-to-flu x ra- 
tio for a sheet-like cloud by Nakano & Nakamural (11978b , the 
critical field strength is given by 



(8) 



with the cloud mass M, and the (projected) area A. The es- 
timates (Table [B are larger by a factor of a few than mag- 
netic field strengths f rom CN Zeeman m easurements (e.g. 
iFalgarone etaD 120081: ICrutcher et all 1201 Oh . suggesting that 
magnetic fields are unlikely to provide wholesale support to 
the IRDCs, although they might be strong enough to affect 
the gas dynamics. 

5. CONCLUSION 

In this paper, we have furthered our analysis of the ammo- 
nia maps presented in Paper 1 , focusing here on IRDC kine- 



matics. Our main conclusions are as follows: 

• In general, the imaged kinematic properties derived 
from the NH3 (1,1) line and (2,2) line are very simi- 
lar and strongly corroborate each other. A notable ex- 
ception is the IRDC G009. 86-0.04 where the line cen- 
ter velocities are offset by ~ 0.5 km s" 1 and the (1,1) 
linewidths are everywhere below 2 km s" 1 while the 
(2,2) linewidths reach up to 4 km s" 1 ; these differences 
are likely due to the active cluster formation underway 
in this IRDC, which selectively affects slightly warmer 
gas traced by the higher excitation (2,2) line. 

• For all of the IRDCs with robust measurements, non- 
thermal motions are greater than thermal motions by 
factors of 2 to 8. The linewidths are always greater 
than the range of centroid velocities across the cloud. 
Indeed, objects of this mass are expected to be in early 
phases of fragmentation from turbulent molecular cloud 
complexes, and this phase of fragmentation is integral 
to setting the conditions of the massive star and cluster 
formation to follow. 

• The velocity fields across the IRDCs are typically very 
regular, showing smooth gradients in centroid velocity 
at the resolved size scales. These gradients could be due 
to rotation, shear, infall, or residual turbulent motions 
from the fragmentation process. Observed departures 
from the regular trends are generally connected to mid- 
infrared point sources tracing embedded young stellar 
objects. At the sites of these sources, the centroid veloc- 
ity may be shifted by 0.5 to 1.5 km s" 1 , perhaps due to 
infall onto, or outflow feedback from, protostars within 
the clouds. These effects tend to be greatest when a 
point source is detected at 24 jim only, i.e. at an early 
phase of star formation. 

• For all of the IRDCs, the kinetic energy estimated 
from the observations is insufficient to provide support 
against collapse. We perform basic models taking into 
account the projected geometry of the IRDC. This spa- 
tial analysis of the thermal, kinetic and gravitational 
energy content indicates that none of the clouds are 
in equilibrium. Rather, the energetics combined with 
the density structure suggest that the clouds are in ac- 
tive fragmentation and collapse, in contr ast to the static 
"turbulent c ore" picture outlined by McKee & Tanl 
(2002, 2003). 
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